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Abstract 

We consider the baryon parity signals at HERA for the case of the MSSM 
production mechanisms and the decays via the lepton number violating couplings 
LiQD. We can probe very small Yukawa couplings A ~ 3 • 10 , limited only by 
the decay length of the LSP. We assume the LSP to be the lightest neutralino and 
study its decays in detail. We present the matrix element squared for the tree- 
level decay amplitude of a generally mixed neutralino explicitly. We find that the 
branching fraction to charged leptons strongly depends on the SUSY parameters 
and can differ significantly from the naively expected 50%. The SUSY mass 
reaches of the studied processes in the ZEUS detector at HERA were found to 
be: (m(e, v)+m{q)) < 170 GeV , 195 GeV and 205 GeV for the L T QD, L^QD and 
L e QD couplings respectively. These are well above existing limits on R-parity 
violating SUSY from previous experiments. We conclude that HERA offers 
a very promising discovery potential for ft p SUSY. 
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1 Introduction 



When extending the Standard Model to include supersymmetry JT| the superpoten- 
tial contains terms which lead to unsuppressed proton decay The most elegant and 
economic solution to this problem is to impose a discrete symmetry which disallows 
the dangerous AB ^ terms. Several discrete symmetries have been proposed, with 
R-parity (R p ) U and Baryon Parity (B p ) [0, Q being the most common. Imposing R p 
results in the minimal supersymmetric standard model (MSSM); in the case of B p the 
superpotential contains the following lepton number violating terms beyond those of 
the MSSM § § 

Kjk-LiQjDk, XijkLiLjEk- (1) 

L, Q denote the left-handed lepton- and quark- doublets, respectively; E, D denote the 
electron and down-like quark singlets respectively. A, A' are Yukawa couplings and 

k = 1,...,3 are generation indices. Theoretically it is still an unresolved question 
which is the preferred or proper symmetry.[] Presumably it will be determined from 
the embedding of the supersymmetric standard model in a more fundamental theory 
at a higher energy scale. In fact, recent work in this direction hints at B p being the 
preferred symmetry 0. 

In light of this theoretical dilemma we propose to consider the full phenomenological 
possibilities, in particular to include the B p signals in all supersymmetric searches. 
There is a growing amount of work in this direction for searches at hadron colliders 



||]-[|ril] and at e + e~-colliders fll2| , [T3| though clearly significantly less than for R p . 

In this paper we focus on promising B p signals at HERA. HERA is particularly 
sensitive to the first operator in Eq.(|I|) since it couples to both leptons and quarks. 
When including a substantial coupling A^- fc the phenomenology changes in two main 
respects: (1) resonant squark production is possible at substantial rates; (2) the lightest 
supersymmetric particle (LSP) is not stable and can decay in the detector. 



Resonant squark production at HERA has been considered in J14 , p~5f| and was 
studied in detail in [TJ| [T7| , where cascade decays of the squarks to the LSP have been 
included. The squark mass reach was found to be fllT 

rriq< 270 GeV for A' > 0.08, (2) 

and the reach in the Yukawa coupling in this mode is 

A'>5.3-10~ 3 for m f ~100GeV. (3) 

In this paper we fully exploit the second possibility, the decay of the LSP. We 
consider the associate pair production of supersymmetric particles via supersymmetric 
gauge couplings 

e~ + q — > e~ + g, (4) 
e~ + q — > v + q' . (5) 



1 Experimentally clearly no statement can be made since supersymmetry has not been found. 
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followed by the decay within the detector of the scalar fermions (/ = e, u, q) via the 
LSP to a .Rp-even final state.Q For example in the case of a dominant operator LiQjD^ 



+ + {(ef + u, + 4) or {yi + dj + 4)}- (6) 

We have denoted the LSP by Xi the conventional notation for the lightest neutralino. 
The main signal for these events will then be a high charged lepton. 

The mechanisms (UH) have been studied in detail in the context of R p , where the 



dominant signals are missing pr [19], |2(J, |21] . These studies have generally lead to the 



conclusion that there is little discovery potential for supersymmetry at HERA beyond 



the limits set by the experiments at the Tevatron. In [|16| |17| this was shown not to be 
the case if Rp is violated. One essential point of this paper is to provide further viable 
signals at HERA for the case of B p . 

In studying the B p signals we go beyond previous work [|T^, [T7|, || in the following 
respects: 

• The production cross sections (f|,[D are independent of the Yukawa coupling. Thus 
it is possible to probe extremely small couplings, which are only limited by the 
requirement that the LSP decay within the detector, i.e. X' ~ 3 • 10~ 6 . Recall, 
the smallest Yukawa coupling in the SM, the electron-Higgs coupling, is just of 
this order. 

• Since the lepton number violating coupling is not involved in the production 
cross section, we are able to explore the full generation structure of the operator 
LiQjDk. We must only require j ^ 3 since a top quark in the final state is 
beyond the kinematic reach of HERA. In particular we can also study the L T QjDk 
operator, which leads to a tau lepton in the final state. 

• We have considered a general neutralino as the LSP (as opposed to restricting the 
study to a photino). This has two main effects in B p studies. (1) For the operator 
LQD the branching fraction for the LSP decay to charged leptons very strongly 
depends on the mixing parameters of the neutralino sector. (2) For a dominantly 
Zino/Higgsino LSP and small LSP masses the decay rate (proportional to a SM 
Yukawa coupling squared) can become so small that the LSP does not decay within 
the detector. This severely limits the range of couplings that can be probed by 
the B p high p? charged lepton signals. 

• We have included the charged current production mechanism ([5]). In Rp this is 
not very promising since it is difficult to distinguish from the SM charged current 
events. However, in B p both the sneutrino and the squark decay via the LSP to 
predominantly visible final states. This increases the overall event rate by about 
a factor 2.5 IE3. 



The clear disadvantage of employing the mechanisms (Q) and fl5|) is that the mass 
reach is severely restricted by the kinematic limit of HERA, since two sfermions must 

2 The direct decay of the sfermions via the operator LQD is strongly suppressed for the small Ay fc 
we are considering. For a detailed discussion see jj], |ll| . 
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be produced. This corresponds directly to the problem encountered in the context of 
the MSSM. However, since we are now considering the case where the LSP decays the 
previous MSSM squark mass bounds derived mainly from searches for missing px signals 
at CDF must be reexamined. We discuss this in more detail below. 
In our analysis we shall make the following simplifying assumptions 

(Al) We shall assume the LSP is a neutralino f\ In some regions of the MSSM param- 
eters the LSP is a chargino. We do not study these supersymmetry models. 

(A2) Of the 27 operators LiQjDk we shall assume in turn that one operator dominates 
[0, and the others are negligible. We thus have 18 separate scenarios to consider 

(A3) The cascade decay of the sfermions to the LSP is 100% of the branching fraction. 

Our paper is outlined as follows. In Section 2 we summarise the relevant previously 
existing bounds on the B p parameters. In Section 3 we review the work on the pro- 
duction mechanisms Eqs.(^J^). In Section 4 we present a detailed analysis of the LSP 
decays for the operator LQD. We focus on the branching fraction to charged leptons 
and on the decay length. Combining the last two sections we examine the signals, 
which we compare with possible backgrounds in Section 5. In Section 6 we formulate 
appropriate cuts to extract the signals. This is the core of the analysis. We discuss the 
results and conclude in Section 7. 



2 Existing Bounds 



The existing bounds on the operators X' i j k LiQjDj C are of two kinds: (a) indirect bounds 
from low-energy processes involving virtual supersymmetric particles; (b) direct bounds 
from collider searches. Both sets of bounds make the simplifying assumption (A2) of a 
single dominant ^-operator. 



(a) The most stringent relevant indirect bounds are |23 
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Note that the bounds scale linearly in the relevant supersymmetric scalar fermion masses 
rh. And there are no indirect bounds on the tau-operators L T QiDj. 

(b) The most detailed collider bounds have been determined by D.P. Roy ||. Here 
use is made of the CDF top-quark search which sets bounds on any di-lepton production 
beyond that of the SM. At the Tevatron, the MSSM pair production of squarks, followed 
by the cascade decay via the LSP would lead to a di-lepton signal. The bound on 

3 This is a strong assumption, since the LSP is not stable and the bounds on relic densities of charged 
or coloured particles are no longer relevant. 
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the di-lepton production directly translates to a bound on the squark pair production 
cross section, which to a good approximation only depends on the squark masses. The 
bound on the squark masses thus obtained is 

m q >100GeV, L e QiDj, L^QiDj. (9) 

And again there is no bound for the tau operators. In addition there are model inde- 



pendent bounds from LEP, which hold for all operators [24] 

m e > 45 GeV, m q > 45 GeV. (10) 

Thus our best bounds are 

m e + m q >iit 5 ^ L -S?' (ID 



90 GeV, L T QD 



Besides studying SUSY bounds at HERA as opposed to the Tevatron, we differ from 
H in the following two important points. First, we consider a general neutralino LSP. 
This has a significant effect on the branching fractions of the LSP in the case of the 
operatoif] LQD. It can also have a significant effect on the lifetime, and thereby the 
decay length of the LSP. The reach in coupling constant in our analysis only depends 
on the decay length (see also Section 4.3). Second, we also consider the tau operator 
L T QD. HERA has an inherently cleaner environment than the Tevatron and thus has a 
non-trivial discovery potential for L T QD. Third, we consider single lepton signals and 
are thus less sensitive to model-dependent variations in the charged lepton branching 
fraction. 

Further collider searches have been proposed for e + e~ |13|] and for hadron colliders 



H, |10|| . In particular also does a detailed analysis of the tau-number violating 
operators for LEP200. They consider gaugino pair production which is complimentary 
to the dominant HERA signal. The signal at LEP200 is insensitive to the squark mass 
and depends on the slepton mass only via the propagator, however it strongly depends 
on the (neutralino) LSP mass. 



3 Production Mechanism 

The dominant R p conserving SUSY production mechanisms at HERA are the Neu- 
tralino Current (NC) and Chargino Current (CC) processes of Figure [I]. Here v, e~L,R 
and qh.R are the sneutrino, the left and right handed selectron and squark respectively. 
In the MSSM there are two pairs of charginos, X121 an d f° ur neutralinos, Xi 4- Follow- 
ing the notation of Haber and Kane [Q], masses and couplings depend on the parameters 
M' and M (the U(l) and SU(2) gaugino parameters), the Higgsino mass mixing param- 
eter \i, and tan(/?) = ^, the ratio of the vacuum expectation values of the two Higgs 
doublets. We shall throughout assume the GUT relation j2~5| 



M' = ^ tan 2 (8 W )M, (12) 



That the operator LQD is special in this respect was already pointed out in |lq, hi 
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where 9\y is the electroweak mixing angle. Gaugino mixing has been thoroughly dis- 
cussed in the literature, and we refer to [[T], [TT| ^IJ for details. 



Processes (|j) and (||) have been calculated by various authors |27|, ^8], |22|| . As 



discussed in [|29||, the cross sections exhibit a strong dependence on the sum of the two 



final state SUSY particle masses, (mg + rriq) or {rrii, + m ? -). Furthermore they show 
some dependence on the gaugino mass parameter M', but are relatively insensitive to 
the parameters [i and tan(/5). Figure |2| shows a contour plot of the SUSY cross section 
& e+q ^i + qi (where I can be a selectron or a sneutrino) as a function of M' and (1111 + rriq). 
The hatched region indicates where the mass of the lightest neutralino becomes larger 
than mi and rriq, the masses of the scalars 5 . Since we have assumed the lightest 
neutralino to be the LSP, we do not consider this region of parameter space. Note that 
within the region of interest to HERA (scalar masses up to ~ 200 GeV) cross section 
variations as a function of M' are within a factor of two for M' > 10 GeV. Figure 
shows the number of expected events per 200 pb^ 1 (two nominal years of HERA running 
for a number of values M' and their corresponding LSP masses. 
The produced (on-shell) selectrons, sneutrinos and squarks subsequently cascade- 
decay into lighter SUSY particles via (making use of assumption (A3) in Section 1) 

el,R + (13) 

y - » + x- (14) 

U,r -> q + xl- (15) 

Eventually only one type of SUSY particle remains, the LSP. In R p conserving models 

the LSP escapes detection, and the resulting SUSY signals are e~+ljet+ missing transverse momentum 

for the NC case, and ljet+pr for the CC case. Several studies |29j have concluded 

that HERA's eventual discovery reach on selectron, sneutrino and squark production 

for R-parity conserving MSSM SUSY is rrig + rriq < 200 GeV and + rriq < 180 GeV 

at an integrated luminosity of 200 pfr -1 . 

In the following sections we shall concentrate on the decay of the LSP via fi p 
couplings, including the full neutralino mixing of the MSSM. We shall show that the 
squark mass discovery reach of the LQD operator is comparable to that of the R-parity 
conserving case, for a very wide range of Yukawa couplings. 



4 LSP Decays 
4.1 Matrix Element 

Here we present the matrix element squared for the decay of the LSP via the operator 
LiQjDk to the charged lepton final state. We consider a general neutralino LSP and we 
retain all final state fermion masses. Thus we also include the "higgsino"-like couplings 
which are proportional to the fermion masses. 

\M(xo^e tUj d k )\ 2 =8c^ 2 Ag fc { (16) 
5 Assuming e, v and q to be degenerate in mass. 

6 Presently, prior to the 1994 running period, the integrated luminosity at HERA is lp6 _1 . 
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D(uj) 2 ei ■ d k [{a{uj) 2 + b(u) 2 )xo • Uj + 2a(uj)b(u)m Uj M% } 
+ D(d k ) 2 ei ■ uj [(a(4) 2 + b{d) 2 )xo • 4 - 2a{d k )b(d)m dk M jio ] 
+ D^jfuj ■ d k [(a(e t ) 2 + b(e) 2 )xo • e; + 2a( ei )6(e)m ei M^ ] 

- D(ei)D(uj) [a(uj)a(ei)m ei m u .Xo • d k + a(uj)b(e)m u .M^ ei ■ d k 

+a(ei)b(u)m ei M^ Uj ■ d k + b(e)b(u) g(uj, xo, e*, d k )} 

- D{uj)D{d k ) [a(uj)a(d k )m u .m dk Xo • - a(uj)b(d)m U:j M^ ei ■ d k 

+a{d k )b(u)m dk M^ei ■ Uj - b(u)b(d) g(uj, Xo, 4, e<)] 

- D{ei)D{d k ) [-a(e i )6(d)m ei M^ M i • d k + a(ei)a(d k )m ei m dk xo ■ uj 

+a(d k )b(e)m dk M< io ei ■ uj - b(e)b(d) g(xo,ei,Uj,d k )]j 

Here Cf = 3 is the colour factor [] and g is the weak coupling constant. We have denoted 
the 4-momenta of the initial and final state particles by their particle symbols. The 
functions D(pi) denote the propagators squared for particle p and are given by 

D(p t )- 1 = Ml + m 2 n - 2x ■ p t - m 2 pi . (17) 



The coupling constants a(pj), b(p) are for example given in |30| and are listed in the 
Appendix for completeness. 

The amplitude squared of the decay to the neutrino, xo —* ^idjd k , can be obtained 
from the above result by a set of simple transformations of the 4-momenta, the propaga- 
tor functions D and the couplings a(pi), b{p) : e, — > z/j, Uj dj. In order to determine 
the total decay rate these two modes must be integrated over phase space, added, and 
then multiplied by a factor of two, since the LSP is a Majorana fermion and can decay 
to the conjugate final states. 

The result for the operators LiLjE k is completely analogous, except the colour factor 
Cf — 1. The result for the operators UDD is given in the Appendix. 



4.2 Charged Lepton Branching Fraction 

Before discussing the results of the calculation in detail, we outline their general nature. 
We focus on the discussion of the LSP decays to r's via the LiQjD k operator (j ^ 3). 
Branching fractions to e ± and ^ (via the operators L e QjD k and L^QjD k respectively) 
are within a few percent of the tau results. The only noticeable difference enters in 
the Higgsino region, where the LSP lifetime is larger for decays to electrons and muons 
than for decays to taus (by a factor {m T /m e ^) 2 ). 

We have found that the branching fraction of the LSP decay to r's and iv's strongly 
depends on the nature (i.e. the gaugino mixing parameters) of the LSP. The LSP decay 
falls into one of the following three distinctively different regions: 

1. A region where the LSP predominantly decays to r's. Branching fractions of the 
LSP in this region are typically BF(x° — > r ± + 2je£s)~70%. Experimentally this 

7 A simpler form of this result was obtained in Ref(l7) for the case of a photino LSP and neglecting 
the final state fermion masses. The colour factor was omitted there by mistake. This error does not 
affect any of the results presented there since only branching fractions were employed, not the absolute 
decay rate. 
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region is the most interesting one, because the event topology contains an "exotic" 
lepton (the r). 

2. A region where the LSP predominantly decays to u T 's. Branching fractions of the 
LSP in this region are typically BF(x1 — > v T + 2jets)~85%. One can still look for 
the signature of decays of the LSP to r's, but the bounds are significantly weaker. 

3. A region in which Higgsino/Zino contributions to the LSP become dominant and 
the LSP mass is small. In this region, the decay rate of the LSP can become so 
small that the LSP only decays within the detector for anomalously large couplings 
A'. Thus the dominant signature would again be the missing j>t signal of the R p 
conserving MSSM scenario. We will not consider this region in the subsequent 
Monte Carlo (MC) study. 

The neutralinos are admixtures of their weak eigenstates 7, Z, and H®. If the 
LSP dominantly consists of the 7 component, branching fractions to are large, while 
a predominance of the Z component makes Xi ~^ v t + 2jets the more preferred decay 
channel. If the LSP is more Higgsino-like, the scenario gets more complicated as the 
branching fractions now also depend on both the vacuum expectation values of the SM 
Higgs and the mass of the decay leptons. 

We shall now discuss the results of the LSP calculation in detail. The amount 
of mixing of 7, Z and the two Higgsinos is determined by the SUSY parameters M', 
M, fi and tan(/3). Assuming the GUT relation (0), we are left with the three "free" 
parameters M', /1 and tan(/3). Figure £| shows the resulting LSP decay regions as a 
function of the SUSY parameters for the dominant operator L T QD. The white region 
indicates where m v o > ±, the region where the mass of the lightest neutralino is 
greater than the mass of the lightest chargino. Since we have assumed that the LSP is 
a neutralino, we do not consider this region in parameter space. The branching fraction 
BF(xi — > + 2jets) is less than 30% | in the light shaded area, and greater than 
30% in the dark shaded regions. We have found that, over a wide range of parameter 
space, the branching fractions change abruptly when the nature of the LSP changes. 
This is demonstrated in Figure [| by solid contours, where the 7, Z and H symbols on 
either side of the contours show which component's contribution dominates the LSP. 
Furthermore we have found that the branching fractions are fairly constant in a given 
branching region, and approximately BF(xi — > + 2jets)~ 70% in the dark regions, 
and BF(xi — > + 2jets)~ 15% in the light regions. 

It is important to notice that these branching fractions significantly modify the 
bounds obtained in Ref . |§ , where as a first approximation a photino LSP was assumed. 



8 We are interested in the pair production of sparticles via processes ([|) and (|[), and the pro- 
duced selectrons, sneutrinos and squarks will subsequently decay to pairs of LSPs (plus other de- 
cay products). The condition that BF(x? — > t ± + 2jets)< 30% corresponds to the condition that 
BF(2x? -> 1 or 2 t^s + jeis)~ 51%. 
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4.3 LSP Lifetime 



In order to get a qualitative understanding of the effects of the neutralino mixing on 
the LSP lifetime let us consider first a photino LSP. A photino LSP decays within the 
detector (defined by the decay length ct^tlsp being less than lm, tlsp is the LSP 
lifetime) provided |7|, |31| 

Kjk > 3 - 7 7 ■ KT 6 (m//100 GeV) 2 (A5 GeV/M LSP ) 5/2 (18) 

where 7 is the Lorentz-boost factor in the laboratory, c is the speed of light. However, 
if the LSP is a pure Higgsino, then the decay rate is proportional to the heaviest mass 
squared normalised to the Higgs vacuum expectation value squared 2m 2 /v 2 , and the 
corresponding bound on A' is 

X' ijk > 8.2 7 • 10- 5 (m r /(m i ))(m / /100G'e\/) 2 (45G'eV r /M L5 p) 5 / 2 . (19) 

For the L^Q\Di operator, the coupling A' ~ 1.4 • 10~ 3 , which would still allow a novel 
region to be explored. For the operator L e QD and a Higgsino LSP, the single squark 
production followed by a direct R p decay should be more promising (discussed in |15|, 
|14|j). If the LSP does decay outside the detector we are left either with the leptoquark- 
like signature or the standard missing pt signature of the MSSM. 

In order to quantify the dependence on the MSSM mixing parameters, we have 
plotted the regions of parameter space (in black) in Figure |] for which the LSP (at a 
value of A' = 0.03) decays outside the detector, for an operator L T QD. These regions 
coincide with regions where the LSP is dominantly a mixture of Higgsinos/Zinos and the 
LSP mass is small. This is shown in Figure [|, which plots the LSP lifetime as a function 
of the Higgsino mass parameter fi. The two peaks at \i ~ GeV and fi ~ 100 GeV 
correspond to the Higgsino/Zino region discussed above. The gap in the middle at 
~ 50 GeV corresponds to the region where the LSP is not a neutralino. Figure |4] also 
shows how far the black regions would extend (dashed lines) at the very small values 
of the coupling constant A' = 3 x 10~ 6 . 



5 Signals and Backgrounds 

In order to investigate the full phenomenological consequences of the LSP decay via the 
LiQjDk operator, we have used a R p violating SUSY generator to simulate the SUSY 
production mechanism of Figures [L]a and [L|b and the subsequent Rp LSP decays at the 
parton level. Throughout our analysis we have chosen to use the MRSD~' structure 



function set of reference |32j for the SUSY events and the background samples. We have 
generated the events at electron and proton beam energies of 26.6 GeV and 820 GeV 
respectively. The SUSY events and potential backgrounds to the SUSY signals were 
passed through a Monte Carlo detector simulation of the ZEUS experiment based on the 
GEANT program [^], and subsequently reconstructed using the ZEUS offline programs. 
We will now proceed to discuss the individual SUSY samples which we have used in 
our analysis, their signatures, and the potential backgrounds to each of the samples. 
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5.1 SUSY Signals 

We have used the MSSM generator [Q to simulate the ftp SUSY samples. In order 



to reduce the large SUSY parameter space, we have made a number of simplifying 
assumptions in addition to assumptions (A1)-(A3): 

• We have assumed sleptons (selectrons and sneutrinos) and the five squarks to be 
degenerate in mass. Throughout we shall denote the mass of the SUSY scalars as 
M-susy = mi = nig. 

• We have only considered one set of gaugino mixing parameters (M' = 40 GeV, 
tan(/5) = 1 and fi = — 200 GeV) in the region where the LSP predominantly 
decays to charged leptons, and one set of parameters (M' = 55 GeV, tan(/3)=4 and 
jj, = +200 GeV) in the region where the LSP predominantly decays to neutrinos. 
As discussed in sections |3] and [|, variations of the gaugino mixing parameters 
within a given LSP decay region have a relatively small effect on the cross section 
and LSP branching fractions (we have found that the combined changes in the 
cross section and the branching fraction are within a factor of two). This is to 
be compared to the large changes in the cross sections as a function of (mj + rriq) 
which are approximately three orders of magnitudes over the SUSY mass scale of 
interest to HERA (cf. Fig. 2). 

We now turn to the signatures of the SUSY events. From the production mechanisms 
of Figure |]a and |I]b two sparticles, a slepton and a squark, are produced. The sleptons 
and squarks subsequently cascade-decay via processes (|T3|)-(^5j) to neutralinos. Let us 
consider two specific examples of events in which the LSP decays via the L T QD operator 
(see also Figure |G|): 

e~ + q -> e~ + q -> e~ + q + 2$ -> e~ + q + r + u T + Ajets. (20) 
e~ + q -> e~ + q -> e~ + q + 2^° -> e~ + q + 2v T + 4jets. (21) 

Processes (ED]) and ( pl|) are both (NC) processes. Note that the pair of sfermions 
will eventually cascade-decay to two LSPs. Process ( ppf) has a tau in its final state. 
This can subsequently decay to a positron or a muon via 

-> e + z/ e z/ T (BF ~ 8.9%) (22) 
-» fi^u^BF ~ 17.8%) (23) 

and the branching fractions of the decays are shown in brackets. Our overall signals 
for LSP decays via the L T QD are thus a positron or a muon and jet activity in the 
detector. For the L e QD and the L^QD operator, our signals are one or more positrons 
or one or more muons respectively and jet activity in the detector]]. 

Having discussed the signals of interest, we now turn to the generated SUSY samples. 
We have generated samples of SUSY events for each of the three operators above. The 



9 Due to the large neutral current deep-inelastic scattering background, the decay of the r to an 
electron is considered a very difficult signature to observe experimentally. 
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± OVJ 


LnQD 


95 


40 


1 


-200 


0.80 


0.16 


400 


32 


L^QD 


100 


40 


1 


-200 


0.80 


0.09 


200 


18 


L T QD 


45 


40 


1 


-200 


0.79 


10.09 


2000 


2018 


L T QD 


65 


40 


1 


-200 


0.79 


2.56 


1000 


512 


L T QD 


80 


40 


1 


-200 


0.79 


0.75 


1000 


150 


L T QD 


85 


40 


1 


-200 


0.79 


0.46 


400 


92 


L T QD 


90 


40 


1 


-200 


0.79 


0.28 


200 


56 



Table 1: Generated R p violating SUSY samples. Msusy is the mass of the scalars; M', 
tan(/?), // are the gaugino mixing parameters; Branch is the branching ratio of the LSP 
to a charged lepton BF(x® — > + 2 jets)] otot is the total cross section of the sample; 
Ngen is the number of events generated; Nexp is the number of events expected for an 
integrated luminosity of 200pb~ 1 , two nominal years of HERA running. 

Le(ti,r)QD sample only contains events in which at least one of the two LSPs has decayed 
to a positron (muon or tau respectively). 

For each type of operator, L e QD, L^QD and L T QD, we have generated MC samples 
at a number of SUSY masses Msusy (= mi — Trig). The samples were generated with the 
gaugino parameters M' = 40 GeV, tan(/?)=l and ji = — 200 GeU, which correspond to 
a region where the LSP predominantly decays to charged leptons. The diagonalisation 
of the gaugino mass matrix produces a LSP with a mass of Mlsp = 44 GeV. Table [I| 
summarises the SUSY samples, their cross sections and the expected number of events 
for an integrated luminosity of 200pb _1 . 

5.2 Backgrounds 

The most promising signatures of the fi p SUSY events are a positron or muon in 
the final state. Backgrounds to such a signature come from "physics backgrounds", 
which produce at least one positron or muon directly, or from "fake backgrounds", in 
which either detector effects fake positron or muon signals, or from positrons or muons 
produced during the fragmentation process from decays of secondary particles. The 
individual background samples are described below. 

In anticipation of a cut in total transverse energy Et in the following analysis, and in 
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Sample Name 


Generator 


Generated Q 2 


range (GeV 2 ) 


(JToripb) 


Ngen 


Nexp 


bb 


HERWIG 


500- 


87248 


4.506 


3395 


901.2 


cc 


HERWIG 


500- 


87248 


37.39 


4105 


7478 


W (e+, m) 


EPVEC 


500- 


87248 


0.219 


400 


43.8 


NC 


LEPTO 


1000- 


+87248 


230.0 


7149 


46000 


CC 


LEPTO 


500- 


-87248 


49.5 


3541 


9900 



Table 2: Generated background samples, otot is the total cross section of the sample; 
Ngen is the number of events generated; Nexp is the number of events expected for an 
integrated luminosity of 200pb _1 . 



order to keep the number of events in our background samples down to a manageable 
size, we have applied a Q 2 > 1000 GeV 2 cut to the deep inelastic scattering neutral 
current sample, and a Q 2 > 500 GeV 2 cut to all other samples at the generator level. 
Kinematically this corresponds to a cut in total Et, Et > 60 GeV, and Et > 45 GeV 
respectively (see Figure ^). However, because a fraction of the transverse energy in each 
event will be missed in the detector and also carried off by neutrinos and hence will not 
be measured, the total visible (or reconstructed) transverse energy will be below the 
kinematically generated Et of Figure |7|. Also calorimeter energy resolution effects will 
broaden the measured Et spectrum, and hence smear events in Et and shift them to 
higher/lower values of Et- Even taking both effects into account, our generator cuts 
on Q 2 are well below our final cut of E T > 90 GeV in the following analysis. 



5.2.1 Physics Backgrounds 

We have considered two types of physics backgrounds. The first one is heavy quark pro- 
duction via high Q 2 photon-gluon fusion [35]. The most dangerous case of photon-gluon 



fusion is expected to be bb production. The produced B° can decay semileptonically 
to a positron or a muon. We have also considered cc backgrounds, but note that these 
are generally less dangerous as their Et distribution is lower. The HERWIG generator 
|36H was used to generate bb and cc events. Table ^| summarises all background sam- 
ples, their cross sections and the expected number of events per luminosity of 200pb~ 1 , 
nominally two years of HERA running. 

The second physics background considered is W production |J7| via the process 

7* + q -> W + q', (24) 

where the virtual photon is radiated from the electron, and the W can decay to give 
a positron or a muon. Cross sections are small compared to the photon-gluon fusion 
backgrounds, but the pt of the charged leptons is in general large. The EPVEC gen- 
erator of reference [[37]] was used to generate the W sample. The sample only contains 
events in which the W decays to a positron or a muon. 
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5.2.2 "Fake Backgrounds" 

Fake signals of positrons or muons in the ZEUS detector can be caused by a number of 
effects, and these include 

• Misidentification of the curvature (and thus charge) of high px electrons in the 
tracking detector. 

• Hadrons depositing energy in the calorimeter with shower profile characteristics 
similar to that of a positron (misidentification of positrons in the calorimeter). 

• Punch-through of high energetic charged hadrons through the calorimeter will 
give rise to tracks in the muon chambers, and hence can lead to misidentified 
muons. 

Furthermore, decays of secondary particles produced during fragmentation processes 
can lead to genuine positrons and muons. In order to investigate all the effects giving 
rise to fake backgrounds, we have generated deep inelastic neutral current (NC) and 
charged current (CC) events using the LEPTO program f38| . Table |] summarises these 
samples. 



6 Analysis 

6.1 Positron and Muon Identification 

Positron and muon identification play a key role in the extraction of the SUSY signals. 
The identification algorithms have to give a high purity in order not to contaminate 
the search sample with too much fake background. With this in mind, we have chosen 
the methods below for positron/muon identification. 

Positron identification is achieved by first of all finding electromagnetic clusters in 
the calorimeter, e.g. by finding shower profiles characteristic of those of electrons. We 
require the cluster to have an energy above 5 GeV. Tracks in the central drift chamber 
(CTD) are then extrapolated to the face of the electromagnetic calorimeter. If a track 
matches such an electron cluster, and the curvature of the track is positive, then the 
cluster is identified as a positron. 

Muon identification is achieved by matching CTD tracks with track segments in 
the muon chambers. Note however that only high energetic muons can penetrate the 
calorimeter before they reach the muon chambers which implies an indirect muon mo- 
mentum cut of p^ > 3 — 4 GeV. 

6.2 Cuts 

Before we go on to discuss how the SUSY signals can be extracted from the background, 
we explain the cuts used in the subsequent analysis: 

• "e + found and /i found" : positron and muon identification have been discussed 
in the previous section. 
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"30 GeV< (E — p z ) < 60 GeV": E is the sum of the total energies of the final 
state particles measured in the calorimeter, and p z is the sum of the z component 
of momentum of all particles measured in the calorimeter. The (E—p z ) variable 
is an important quantity which characterises the event PDA: for e-p events where 



all final state particles are measured, (E—p z ) = 2E e , where E e is the energy of the 
initial state electron. Undetected particles which are emitted down the forward 
beampipe give a negligible loss in (E—p z ), while for example photoproduction (i.e. 
7-p scattering) processes in which the scattered electron remains in the beampipe, 
predominantly result in low values of (E —p z ). The (E —p z ) cut is very efficient 
in rejecting photoproduction events and other backgrounds to the deep inelastic 
scattering processes of interest. 

In our study we have used the (E — p z ) cut to mainly reject charged current (CC) 
events and W events. In both types of background a considerable fraction of the 
energy is carried off by neutrinos, and hence will not be detected. As a result the 
(E — p z ) distribution of such events is shifted towards lower values. 

• "circularity > 1": One can define a 2- dimensional- like sphericity, circularity 
|K| , which is defined in the range <circularity< 1. The circularity (transverse 
sphericity) is defined by 

circularity = 2\ 2 (25) 
in terms of the smaller eigenvalue, A2, of the two dimensional sphericity tensor 

S a /3 = PiaVip/ Pi 2 ( 26 ) 
i i 

where a, (5 are Cartesian components, i runs over the number of particles (i.e. 
reconstructed calorimeter clusters) in the event, and p are the four-momenta of 
the particles. The circularity reflects the isotropy of the event. A circularity of 
unity corresponds to a completely isotropic event. A three jet event, for example, 
with the three jets pointing in directions of the azimuthal angle, <ft, <ft = 0°, 
cf) = 120° and <p — 240° would be a good example of an event with a high value of 
circularity. A two jet event, with a back-to-back jet configuration has a circularity 
close to zero. We shall use this cut to distinguish between the SUSY events and 
backgrounds. As will be discussed below, the SUSY events are very isotropic, and 
thus have a high value of circularity. 

• u Et > 90 GeV" : E? is the total transverse energy measured in the calorimeter. 
We also include the measurement of the muon track in the E? if a muon was 
found. 

6.3 Signal Extraction 

We first examine the characteristics of the positrons and muons in the events. Figure 
|8|a and Rib show the muon and positron px spectra for the SUSY samples at Msusy = 
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The z-axis is defined to point along the proton direction. 
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80 GeV, correctly normalised with respect to each other. Note that the muon and 
positron p T spectrum of the L T QD sample is much softer than the L^QD and L e QD 
Pt spectra. This is because the r's have to decay via processes fl22|) and (|23|) to positrons 
and muons, and hence the energy of the e + or \i is shared with the associated neutrinos. 

Even though the pt spectrum of the L^QD sample is slightly harder than the pt 
spectrum of the backgrounds of section |572|, we conclude that the positron or muon 
Pt variable is not a good discriminator between the SUSY signals and the background 
processes. 

Next we consider the circularity of the events. The decay process ( POD gives a typical 
SUSY event topology, consisting of 5 jets, a r lepton, an electron, neutrinos and the 
proton remnant. Unfortunately not all of the jets are well resolved - most of the jets 
will generally be very soft. We have found that jet finding algorithms can wrongly 
recombine the softer jets with harder jets, and the jet multiplicity of such events can 
be incorrectly determined. We have run jet-finders on the MC samples, and have 
found that the jet multiplicity can be used to discriminate between the SUSY events 
and backgrounds. However, we have also found that the circularity variable is vastly 
superior to jet multiplicity. This can be explained as follows: 

• The high number of jets result in a more isotropic event structure for the SUSY 
events. This makes the events more circular than their respective backgrounds. 

• Soft jets can not be resolved by jet-finders, but they do contribute to the circularity 
of the event. 

Figure |] shows the circularity of all combined backgrounds and the SUSY L T QD 
samples with M S usy = 45 GeV and M S usy = 80 GeV. We observe two features: 

• The circularity of the background falls off much more rapidly than the circularity 
of the SUSY events, therefore a cut in circularity (chosen to be at 0.1) will be a 
good discriminator between the SUSY signal and background. 

• Higher SUSY scalar masses have the effect of making the events more circular - 
as a result higher Msusy samples will be more efficient in passing the circularity 
cut. 



Figure [T0| a shows the E T distribution of the backgrounds with and without a circu- 
larity cut. The cut reduces the background by a factor of 10 at Et = 60 GeV, and by 
a factor of 18 at Et = 140 GeV. A similar plot is shown for the SUSY L T QD samples 
with SUSY masses M S usy = 45 GeV and M S usy = 80 GeV. Again note that for the 
higher mass sample the efficiency of the number of events surviving the circularity cut 
at large values of Et is higher than for the lower mass sample. 

After requiring an e + or /i and applying (E — p z ) and circularity cuts to the MC 
background samples, we are left with (102.4 ± 20.9) events. Figure [TTl a shows their Et 
distribution. The right-hand side of the Et tail is dominated by bb and cc events. Note 
that due to the higher cross section the cc sample is predominant in the Et background 
distribution up to an E T of 80 GeV. The background events with E T > 80 GeV are 
all bb events, which generally have a much harder Et spectrum. We have chosen our 
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final cut at Et > 90 GeV. This includes (0.53 ± 0.31) bb events. Figure |ll| b shows for 
comparison the E T distributions of the three SUSY samples L e QD, L^QD and L T QD 
at high mass ranges. 



6.4 Discovery Reaches 

We now apply the analysis described above to all MC samples. The only difference 
in the cut sequence applied to the L e QD, the L^QD and the L T QD samples is the 
positron/muon selection. We demand 

• One or more found positrons for the L e QD sample. 

• One or more found muons for the L^QD sample. 

• One or more found positrons or muons for the L T QD sample. 

Tables [| f| and || show the effects of the cuts on the SUSY samples and the back- 
grounds. Figure Oa shows the efficiencies of the samples after cuts. Normalising the 



cross sections of table |T| to an integrated luminosity of L=200pb _1 , and folding in the 
efficiencies of Figure |T^ a we obtain Figure T5b. We then use Poisson statistics to in- 
vestigate the potential discovery limit of the L e QD operator at the 90% confidence 
level, by requiring at least 2.3 signal events and assuming zero background. The above 
condition holds over the SUSY mass range 90 GeV < (rrij + rriq) < 205 GeV for the 
L e QD operator and the chosen set of gaugino parameters (M' = 40 GeV, tan(/3)=l 
and = —200 GeV) . The lower mass point at (mj + rriq) = 90 GeV corresponds to the 
current model-independent LEP bounds. Using the same argument for the L^QD and 
the L T QD sample, but requiring at least 3.1 signal events and assuming 0.53 background 
events we can determine the approximate 90%CL discovery reaches for the L^QD and 
the L T QD operators. Note that the gaugino parameters of the above SUSY MC samples 
were chosen in a region where the LSP predominantly decays to charged leptons. 

We have also investigated a region in SUSY parameter space where the LSP predom- 
inantly decays to neutrinos. Table ^ lists the SUSY cross sections and LSP branching ra- 
tios in this region. Because the parameters M' = 55 GeV, tan(/3)=4 and \i = +200 GeV 
produce a LSP with M LS p = 44 GeV, the same LSP mass as the generated MC samples 
of table 0, we can use the efficiencies of Figure 12 a to obtain the expected number of 
events after cuts of Figure |T2|c. 

Using the criteria previously described we now summarise the maximum discovery 
reaches for the three operators L e QD, L^QD and L T QD in the two regions of gaugino 
parameter space: 
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< 


205 GeV for dominant L e QD, 


LSP decays to e ± , 
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m 


% 


y) 


+ m(q) 


< 


185 GeV for dominant L e QD, 


LSP decays to v e ±, 


(28) 


m 
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< 


195 GeV for dominant L^QD, 


LSP decays to [x^, 
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m 


% 


z) 


+ m{q) 


< 


175 GeV for dominant L^QD, 


LSP decays to i/„±, 


(30) 


m 
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+ m(q) 


< 


170 GeV for dominant L T QD, 


LSP decays to r ± , 


(31) 


m 


% 


y) 


+ m(q) 


< 


140 GeV for dominant L r QD, 


LSP decays to v T ±. 


(32) 
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cut 


bb 


cc 


w 


NC 


cc 


none 


3395 


4105 


400 


7149 


3541 


e + found 


22 


24 


118 


55 


8 


30 GeV < (E- Pz ) < 60 GeF 


20 


21 


27 


54 


2 


circularity > 0.1 


4 


7 


9 


5 





E T > 90 GeV 

















Exp. no. of evts 


















cut 


L e QD^ 


L e QDQ5 




L e QD^ 


L e QDi05 


none 


940 


500 


400 


100 


100 


e + found 


423 


230 


177 


38 


44 


30 GeV < (E — p z ) < 60 GeV 


385 


201 


151 


35 


36 


circularity > 0.1 


255 


174 


136 


33 


33 


E T > 90 GeV 


67 


124 


124 


32 


32 


Exp. no. of evts 


95 ± 10 


84.5 ± 6.5 


31.0 ± 2.3 


7.0 ± 1.0 


1.9 ± 0.3 



Table 3: The cut-sequence applied to the MC samples. Here u L e QD^ n refers to the 
SUSY sample with ~M.su sy = 45 GeV, and LSP decays via the L e QD operator. The 
last row shows the expected number of events, properly normalised, after cuts for an 
integrated luminosity of 200 pb~ x . 



cut 


bb 


cc 


w 


NC 


CC 


none 


3395 


4105 


400 


7149 


3541 


[X found 


126 


71 


65 


29 


11 


30 GeU < (E- Pz ) < 60 GeV 


114 


59 


24 


20 


3 


circularity > 0.1 


42 


18 


3 


2 





E T > 90 GeV 


2 














Exp. no. of evts 


0.53 ± 0.31 















cut 


L^QD 45 










none 


1000 


600 


400 


400 


200 


/j, found 


395 


191 


114 


119 


50 


30 GeU < (E- Pz ) < 60 GeV 


330 


157 


92 


102 


44 


circularity > 0.1 


177 


131 


78 


88 


31 


E T > 90 GeV 


27 


71 


58 


79 


31 


Exp. no. of evts 


54 ± 10 


60.7 ± 6.7 


21.8 ± 2.7 


6.3 ± 0.6 


2.8 ± 0.5 



Table 4: The cut-sequence applied to the MC samples. Here a L^QD A ^ n refers to the 
SUSY sample with M S usy = 45 GeV, and LSP decays via the L^QD operator. 
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cut 


bb 


cc 


W 


NC 


cc 


none 


3395 


4105 


400 


7149 


3541 


e + or n found 


148 


95 


183 


84 


19 


30 Gey < (E- Pz ) < 60 GeV 


134 


80 


51 


74 


5 


circularity > 0.1 


46 


25 


12 


7 





E T > 90 GeV 


2 














Exp. no. of evts 


0.53 ± 0.31 















cut 


L T QD 45 


L T QD 65 


LrQD 8 o 


L T QD 85 


L t QDq 


none 


2000 


1000 


1000 


400 


200 


\x found 


185 


74 


92 


33 


14 


30 Gey < (E — p z ) < 60 GeV 


166 


53 


71 


27 


11 


circularity > 0.1 


105 


40 


58 


21 


8 


E T > 90 GeV 


17 


26 


44 


17 


8 


Exp. no. of evts 


17.2 ± 4.2 


13.4 ± 2.6 


6.6 ± 0.9 


3.9 ± 0.9 


2.2 ± 0.8 



Table 5: The cut-sequence applied to the MC samples. Here u L T QD4^ n refers to the 
SUSY sample with M S usy = 45 GeV, and LSP decays via the L T QD operator. 







SUSY 


parameters 








Operator 


^SUSY 


M' 


tan(/3) 


A* 


Branch 


otot 




(GeV) 


(GeV) 




(GeV) 


(%) 


(Pb) 


L e QD 


45 


55 


4 


+200 


0.16 


1.67 


L e QD 


65 


55 


4 


+200 


0.16 


0.43 


L e QD 


80 


55 


4 


+200 


0.16 


0.12 


L e QD 


95 


55 


4 


+200 


0.16 


0.03 


L e QD 


105 


55 


4 


+200 


0.16 


0.01 


Lfj,QD 


45 


55 


4 


+200 


0.16 


3.12 


L^QD 


65 


55 


4 


+200 


0.16 


0.81 


L^QD 


80 


55 


4 


+200 


0.16 


0.24 


L^QD 


95 


55 


4 


+200 


0.16 


0.05 


L^QD 


100 


55 


4 


+200 


0.16 


0.03 


L T QD 


45 


55 


4 


+200 


0.16 


3.12 


L T QD 


65 


55 


4 


+200 


0.16 


0.81 


L T QD 


80 


55 


4 


+200 


0.16 


0.24 


L T QD 


85 


55 


4 


+200 


0.16 


0.15 


L T QD 


90 


55 


4 


+200 


0.16 


0.09 



Table 6: Generated R p violating SUSY samples in a region of gaugino parameter space 
where the LSP predominantly decays to neutrinos. 
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Note that the discovery reaches d27|)-(p2D only hold for the assumptions we have 
made in section |5.1| , and only for the particular choice of the gaugino mixing parameters, 
namely (M f = 40 GeV, tan(/?)=l and /x = —200GeV) for dominant LSP decays to l^, 
and (M' = 55 GeV , tan(/3)=4 and \i = +200 GeV) for dominant LSP decays to v l± . We 
have also investigated the effect of smaller LSP masses on the detection efficiencies of 
the samples, but found that they are generally small: for a LSP mass of Mlsp = 20 GeV 
the correction to the efficiencies of Figure |T2| a (at Mlsp — 44GeV) were found to be 
approximately greater than 0.6. Furthermore, because cross section variations as a 
function of M', tan(/3) and \x are relatively small (see section |3[) and changes in the 
LSP branching fractions to charged leptons are fairly constant in a given characteristic 



region of the LSP (see section |j) we conclude that the discovery reaches of ([27|)-(|32]) 
are approximately valid over a wide range of gaugino mixing parameter space within a 
given characteristic region of the LSP. 



7 Discussion and Conclusions 

We have investigated decays of the LSP as a general neutralino via the LQD operator. 
We have found that the results depend on the mixing parameters of the neutralinos. 
Three distinct regions in the parameter space characterise the LSP decay: (1) A region 
in which the photino component of the LSP is dominant. In this region branching 
fractions to charged leptons are high, typically ~ 70%. (2) A region in which the 
zino component of the LSP is dominant. In this region branching fractions to charged 
leptons are low, and typically ~ 15%. (3) A region in which the LSP is dominantly a 
Higgsino/Zino admixture at low LSP masses. In this region the decay rate of the LSP 
can become so small that the LSP only decays within the detector for anomalously large 
couplings A', in which case the only distinguishing feature of is the single sfermion 
production. In all other cases the LSP decays in the detector for Yukawa couplings 
as low as 10~ 3 . At coulings A' = 3 ■ 10~ 6 the LSP decays inside the detector in only half 
the mixing scenarios. 

Based on the above results we have investigated the discovery reaches of HERA 
on the LQD operator using a full MC simulation of pair produced SUSY signals with 
subsequent $ p decays of the LSP. We have considered the regions (1) and (2) of Section 
4.2 in gaugino parameter space, and have found that fi p SUSY can be discovered at 
HERA up to combined SUSY masses (m(e, 9) + m(q)) of 

• 170 GeV, 195 GeV, 205 GeV for L T QD, L^QD, L e QD and mixing scenario (1) 

• UOGeV, 175 GeV, 185 GeV for L T QD, L^QD, L e QD and mixing scenario (2) 

after two years of HERA running at nominal luminosity. The discovery reaches are to 
be compared with existing limits of 90 GeV, 145 GeV, 145 GeV (L T QD, L^QD, L e QD). 
The limits will only hold for the investigated LSP mixing scenarios and for the set of 
assumptions we have made on the SUSY model, namely that (a) the LSP is a neutralino, 
that (b) only one of the 27 operators LiQjD k is dominant, that (c) sleptons and squarks 
are degenerate in mass, and that (d) the branching fraction of sfermion cascade decays 
to the LSP is 100%. 
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Because HERA is an inherently cleaner environment than the Tevatron, there are 
clear advantages for the search of /2 P SUSY signals. Firstly our analysis demands 
only one charged lepton from one of the LSP decays. Consequently future limits from 
HERA will be less model dependent, and thus complimentary to existing limits from 
the Tevatron 0. Secondly our analysis is sensitive to LSP decays to r's. HERA is the 
first collider experiment which can directly probe the L T QD. The discovery potential 
of r signals at HERA nicely ties in with existing model independent limits on sfermion 
masses from LEP and proposed searches for r signals at LEP 200. 

We conclude that HERA offers a very promising discovery potential for direct 
searches for fi v SUSY, which is nearly twice as high as the existing indirect bounds 
on the L r QD operator. 
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Appendix 

Here we complete the discussion on the LSP decay matrix element. The relevant cou- 
pling constants are given in Table [7|. Here a(pi) denotes the coupling to the Higgsino 
which is proportional to the mass, and b(pi) is the coupling to the gauginos. N(l), N(2) 
are the LSP admixtures of the Bino and the neutral Wino respectively. iV(3), iV(4) are 
the Higgsino admixtures. 

The matrix element squared for decay via the baryon number violating operator 
(which we do not make use of in this paper) is given by. 

\M(xo ^ udA)\ 2 = 8c f g 2 \" 2 { (33) 
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The notation is as in Section 4.1. The colour factor cy = 6. 

11 The analysis of reference g assumes dilepton signals coming from the ft p decay of two LSPs. 
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a(Pi) 


KPi) 


e» 

Mi 

di 
di 


^ 771 

2M W cos /3 " Le i 

N(3) m 

2M W cos & " Le i 



JV(4) 
2M w s\n(i mu ] 

"^(^) 777 

2M w sin/3" tu i 

2M w cos(3" i(l k 

JV(3) 
2M W cos I3 md k 


-±(N(2)+tan9 w N(l)) 

-tan ^yJV(l) 
|(JV(2) -tan^JV(l)) 
±(iV(2) + l./3.tan0 w iV(l)) 

§tan0 w iV(l) 
-|(iV(2) - §tan0 w 7V(l)) 
-|tan%7V(l) 



Table 7: The Coupling constants a(pi) and used in the LSP decay calculation. 
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